Effects of interaction on the diffusion of atomic matter waves 
in one-dimensional quasi-periodic potentials 
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We study the behaviour of an ultracold atomic gas of bosons in a bichromatic lattice, where 
the weaker lattice is used as a source of disorder. We numerically solve a discretized mean-held 
equation, which generalizes the one-dimensional Aubry- Andre model for particles in a quasi-periodic 
potential by including the interaction between atoms. We compare the results for commensurate 
and incommensurate lattices. We investigate the role of the initial shape of the wavepacket as 
well as the interplay between two competing effects of the interaction, namely self-trapping and 
delocalization. Our calculations show that, if the condensate initially occupies a single lattice site, 
the dynamics of the interacting gas is dominated by self-trapping in a wide range of parameters, even 
for weak interaction. Conversely, if the diffusion starts from a Gaussian wavepacket, self-trapping is 
significantly suppressed and the destruction of localization by interaction is more easily observable. 

PACS numbers: 03.75.Lm, 03.75.Kk 



I. INTRODUCTION 

Localization induced by disorder has been recently ob- 
served in_^ultracold bosonic gases in purely random po- 
tentials [1] and in bichromatic optical lattices 0. In the 
experiment of Rcf . [Ij , a Bose-Einstein condensate of in- 
teracting atoms is released from a trap into a speckle 
potential and the localization effects appear in the low 
density tails of the spatial distribution of atoms, where 
the interaction is negligible. In [2], the interaction be- 
tween atoms is suppressed from the very beginning by 
tuning the s-wave scattering length to zero by means of 
a Feshbach resonance. In both cases, the observations 
have been interpreted in terms of Anderson localization 

i- 

An interesting feature of one-dimensional (ID) bichro- 
matic lattices is that they allow one to implement the 
Aubry- Andre model [3| , also known as Harper model Q , 
for particles in a quasi-periodic potential. In the case 
of noninteracting particles, this model predicts a sharp 
transition from diffusion to localization for a given value 
of the disorder strength. In ultracold gases, thanks to 
the availability of Feshbach resonances, the interaction 
between atoms can be changed almost at will, thus allow- 
ing for the investigation of the role played by interaction 
in the transition from diffusion to localization 6]. The 
effects of interaction have been recently the subject of 
several theoretical investigations in the case of localiza- 
tion in purely random potentials J,, &, U^, [HI, [13] and 
quasi-periodic potentials [l^ fl^-llSj. but some results 
are still controversial. It is worth mentioning that the 
Aubry- Andre model has been recently implemented also 
in experiments with diffusion of light in photonic lattices 

The purpose of this paper is to understand the inter- 
play between diffusion and localization of an interact- 
ing Bose-Einstein condensate in a ID bichromatic lattice. 



One of the main problems arises from the occurrence of 
two competing effects of the interaction: on the one hand 
it favours localization through the self-trapping mecha- 
nism fr?! , on the other hand it is expected to destroy the 
(Anderson) localization induced by disorder. In order 
to shed light on these effects and make a bridge between 
theoretical results and feasible experiments, we study the 
following issues: i) what happens when the two lattices, 
which are superimposed to create the bichromatic lat- 
tice, have commensurate wavelengths (periodic system) 
instead of incommensurate (quasi-periodic system, as in 
the Aubry- Andre model); ii) what is the role of the ini- 
tial shape of the condensate; iii) what is the role of self- 
trapping processes. All these issues have subtle implica- 
tions in the experimental observability of the crossover 
between diffusion and localization in an interacting gas. 
Our calculations show that, if the condensate initially oc- 
cupies a single lattice site, the dynamics of the gas in the 
lattice is dominated by self-trapping in a wide range of 
parameters, even for weak interaction. Conversely, if the 
diffusion starts from a condensate with Gaussian shape, 
extended over several lattice sites, self-trapping is signif- 
icantly suppressed and the destruction of localization by 
interaction is more easily observable. 



II. THE MODEL 

One-dimensional bichromatic lattices are realized in 
experiments with Bose-Einstein condensates by superim- 
posing two optical lattices of different wavelengths, pro- 
ducing an external potential acting on the atoms in this 
form: 

V{x) = siEji^ cos^ikix) + S2Eji^ cos^{k2x) , (1) 

where Er- = fi^A:|/(2m) is the recoil energy and Sj is the 
dimensionless lattice strength. One of the two lattices is 
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used as the main periodic potential (primary lattice) fix- 
ing the Bloch band structure of the single-particle states 
without disorder. It is usually strong enough (si 3> 1) 
to apply the tight-binding approximation, i.e., the atoms 
occupy the sites of the primary lattice and can tunnel 
from one site to the other with a given tunneling rate 
J [H, [l^. The second lattice is weaker (s2 <C si) and 
introduces a "deterministic" disorder, or quasi-disorder 
[3, [13, . For a noninteracting gas in such a potential, 
the evolution of the system is described by the Aubry- 
Andre model [3] , which is obtained from the Schrodinger 
equation by expanding the single-particle wavefunction 
tjj{x) over a set of Wannier states, maximally localized at 
the minima of the primary lattice in the lowest Bloch 
band, \^p) = ^jipj\wj) [H, US]- In the presence of 
interactions between the atoms, one can instead start 
from the Gross-Pitaevskii (GP) equation [H, [l^l and use 
the same procedure in order to get a generalized Aubry- 
Andre model which includes an additional nonlinear term 
that represents the mean-field interaction. The Hamilto- 
nian is 



H 



with 



i 



Vj = X coa{2TTaj 



(2) 



(3) 



where j is the primary lattice site index, a = k2/ki is 
the ratio between the wavenumbers of the two lattices, 
9 is an arbitrary phase, Tpj is a complex variable whose 
modulus square gives the probability of finding a particle 
at the lattice site j. In deriving the Hamiltonian ^ from 
the GP equation for a realistic three-dimensional conden- 
sate, we assume the condensate to be axially symmetric 
and elongated in the direction of the lattice, such that 
the transverse confinement is sufficiently strong to freeze 
out the radial dynamics and the longitudinal confinement 
sufficiently weak to be neglected. We have also chosen 
En-^J = 1. The dimensionless parameters A and (3 rep- 
resent the strength of the disorder and of the mean-field 
interaction, respectively, and are the key parameters of 
the present work. 

The equations of motion are generated by id^jjj/dt = 
dH/dil)*, yielding 



^,_i+y,V,+/3|^,fV', (4) 



that can be considered as the GP equation on a dis- 
cretized lattice. Similar versions of a discretized GP 
equation have been already used to investigate the dy- 
namics of condensates in different situations (see for in- 
stance Ref. [13] )• In the above equation, the time t is 
expressed in dimensionless units. The actual time in sec- 
onds can be obtained by multiplying t by h/{JEfi-^). 

We study the localization properties of the system by 
considering the problem of quantum diffusion of an ini- 
tially localized wavepacket. We solve Eq. ^ by using 



a standard fourth order Runge-Kutta algorithm for the 
numerical integration. The accuracy of the integration is 
checked by monitoring the conservation of the norm of 
the wavepacket and of the energy of the system. We in- 
vestigate the evolution starting from two different classes 
of initial conditions, namely a (5-function localized in a 
single lattice site, 



and a Gaussian wavepacket of width cr. 



(5) 



(6) 



where TV is a normalization factor that has to be deter- 
mined according to IV'iP = 1- Owing to arbitrariness 
of the phase 9, here we have chosen, without any loss of 
generality, to fix the initial localization center at j — 0. 

As a measure of the localization we consider two quan- 
tities: the width of the wavepacket measured as the 
square root of the second moment of the spatial distri- 
bution \ipj{t)\'^, 



and the participation number 



Pit) 



(8) 



which measures the number of significantly occupied 
lattice sites [1^. The quantity (j) represents the av- 
erage over the spatial distribution |-!/)j(i)p, defined as 



III. NONINTERACTING PARTICLES 

When /3 = (noninteracting gas) and a is irrational 
(quasi-periodic lattice), Hamiltonian ^ coincides with 
the well-known Aubry- Andre model that has been widely 
studied in the literature (see [18|] and references therein). 
Aubry and Andre [3] analytically proved that the system 
undergoes a transition from extended to localized states 
at A = 2. All states are extended for A < 2, critical for 
A = 2 and exponentially localized for A > 2, with the 
same localization length l/ln(A/2). Aubry and Andre 
A] identified this exponential localization as the Ander- 
son localization in a quasi-periodic potential, analog to 
Anderson localization in a purely random potential [26| . 

From the numerical viewpoint, the static spectrum of 
the Aubry- Andre Hamiltonian can be calculated by solv- 
ing the following eigenvalue problem 

— V-'j+i - V-'j-i + Acos(27raj + 9)ijjj ~ Eipj (9) 

and considering a sequence of rational numbers a„, that 
converges to an irrational a as n — > oo (see for instance 
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[27l.[2^). The sequence of approximants a„ can be found 
by successive truncations of the continued-fraction ex- 
pansion of a. A standard choice consists of choosing the 
inverse golden mean a — (-\/5 — l)/2 25]. In this case 
the approximants are obtained by writing a„ = pn/qm 
where p„ and g„ ~ Pn+i are two consecutive terms of the 
Fibonacci sequence {pi = P2 = ^, Pn = Pn-i + Pn-2 for 
n > 2). The incommensurate case can thus be considered 
as the hmit of a sequence of commensurate Hamiltonians, 
whose eigenvalues E'^'™ and eigenfunctions cj)^'™' can be 
labelled by the quasi-momentum k and the band index 
TO, since the spatial periodicity of the system, with pe- 
riod Qn, permits to use the Bloch wave decomposition. 
One finds that, for sufficiently large n and for A > 2, 
the eigenfunctions are indeed characterized by periodic 
replica of exponentially localized functions within each 
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FIG. 1: Diffusion of a noninteracting gas (/3 = 0) in the 
Aubry- Andre model with a — (\/5 — 1)/2. The time evolution 
of the width of the wavepacket is shown for different val- 
ues of the disorder strength, A = 1.5, 1.7, 1.9, 2, 2.1, 3. In the 
top panel, the starting wavepacket is a (5-function localized in 
a single site. In the bottom panel we use an initial Gaussian 
wavepacket of width o" = 5. In both cases, one clearly observes 
the transition from extended to localized states that occurs 
at A = 2. Here and in the following figures, time and width 
are expressed in the dimensionless units defined in section Ull 
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FIG. 2: (color online) Time evolution of the width of the 
wavepacket wii) of noninteracting particles, starting from a 
single-site ^-function, for A = 2 and for different orders, n, 
of the approximants in the Fibonacci sequence. The black 
arrows represent the values of t at which we observe the tran- 
sition from the behaviour predicted for a quasi-periodic po- 
tential (incommensurate lattice) to the diffusion expected in 
a periodic potential. 

period of the potential, that in the limit n — > c» tend to 
a single localized function 18]. 

The localization transition at A = 2 can be observed 
also in the dynamics (quantum diffusion), by looking for 
example at the width of the wavepacket as a function of 
the time [2§]. In particular, in the quasi-periodic lattice, 
the asymptotic spreading of the wavepacket width w(t) 
can be parametrized as w(i) ~ i''', and one finds three 
different regimes: 

(i) A < 2: ballistic regime, 7 = 1 

(ii) A = 2: sub-diffusive regime, 7 ^ 0.5 

(iii) A > 2: localized regime, 7 = . 

Our results for a — {\fh — l)/2 are shown in Fig. [TJ In 
the case of an initial (5-function wavepacket (top panel), 
we find perfect agreement with previous calculations 29| . 
The lower panel shows our results for the case of an initial 
Gaussian wavepacket. By comparing the two cases, one 
can see that the asymptotic behaviour is not affected by 
the choice of the initial shape of the wavepacket. 

It is worth stressing that a truly quasi-periodic poten- 
tial can not be realized in any realistic experiment with 
ultracold atoms, since the system is finite and the ratio 
of the frequencies of the two laser beams is a rational 
number. It is thus important to clarify to which extent 
the predictions of the Aubry- Andre model are relevant 
for the description of current experiments in bichromatic 
lattices. To this purpose we compare the results of the 
quasi-periodic potential with those obtained by using ap- 
proximants of order n. For the value a = (\/5— 1)/2, this 
consists of selecting the term of order n in the Fibonacci 
sequence. For any finite value of n the system is periodic, 
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FIG. 3: (color online) Modulus square of the wavefunction 
jf/'jp for different values of n, plotted at a fixed evolution 
time t = fOOO, for A = 7 and (3 = 0. The initial wavepacket 
at t = is a (5-function localized at j = 0. The vertical arrows 
are drawn at the positions q„/2. 



with wavelength g„, and the diffusion of an initially lo- 
calized wavepacket is expected to be ballistic {w{t) ~ t). 
However, in the limit n — > cxd one must recover the re- 
sults of the Aubry- Andre model, with a critical behaviour 
for A = 2 and localized states for A > 2. The approach 
to this limit in nontrivial and involves the characteristic 
time and length scales of the system. 

In Fig. [2] we first show our results for the diffusion of a 
(5-like wavepacket in a lattice with the critical value A = 2. 
For any finite n the wavepacket exhibits a sub-diffusive 
spreading (wit) ~ with 7 « 0.5), as in the incommen- 
surate case, within an initial time interval. Then, at time 
r, the width starts growing as in a ballistic diffusion in a 
periodic lattice. The transition between the two regimes 
turns out to occur when the width of the wavepacket be- 
comes of the same order of the spatial periodicity of the 
lattice. The transition time, r, indicated by the arrows 
in Fig. [21 increases with the order n of the approximants 
and the corresponding width, wIt) exhibits a linear de- 
pendence on the periodicity of the system, q„ (s^j . 

The role of the spatial periodicity is even more evident 
if one plots the density distribution in the regime of lo- 
calization, as shown in Fig. O for A = 7 and t — 1000. In 
this figure the arrows are drawn at the positions qri/2. As 
one can see, the deviations from the density distribution 
of the incommensurate case (n — > 00) are caused by the 
spreading of the lateral components of the distribution, 
i.e., those at a distance of the order of, or larger than 
q„/2. The asymptotic behaviour (t — > 00) is always bal- 
listic. However, for a finite t and for A > 2 the central 
part of the density distribution (within a width of order 



Qn) exhibits an exponential localization, independent of 
n, and is almost indistinguishable from the one predicted 
for the incommensurate lattice. The spreading of the low 
density tails affects the behaviour of the width defined in 
Eq. (O. An example is shown in Fig.[4l For short times 
the contribution of the expanding tails is negligible, while 
for later times the width increases as in a ballistic diffu- 
sion. It is worth stressing, however, that these effects of 
the low density tails are expected to be hardly detectable 
in actual experiments, due to the finite resolution in the 
measurement of the density distribution. 

Given the typical timescale [Slj and optical resolution 
of the experiments with ultracold gases in optical lattices, 
our analysis confirms that the transition from diffusion 
to localization observed in Ref. [3] can correctly be in- 
terpreted in terms of the predictions of the Aubry- Andre 
model. 



IV. INTERACTING PARTICLES 

In this section we consider the interacting case (/3 ^ 0). 
We mainly focus on two effects of the interaction, namely 
the self-trapping phenomenon and the delocalization in- 
duced by the interaction in the regime A > 2. As al- 
ready mentioned in the introduction, these two effects 
are competing and must be carefully analysed in order 
to correctly interpret the expansion of a wavepacket. 



A. Self-trapping 

Self-trapping is a localization phenomenon, different 
from Anderson localization, that occurs when the inter- 
action is stronger than a critical value Be, even for a 
purely periodic system without disorder [H, [H, [s^l 
and double- well potentials [H, [s^ . An intuitive under- 
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FIG. 4: (color online) Time evolution of the width of the 
wavepacket w{t) of noninteracting particles, starting from a 
single-site (5-function, for A = 3 and for different values of n. 
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FIG. 5: (color online) Width w{t\ participation number P(t), 
and density distribution \il}j(t = 1000)^ for two values of the 
interaction strength fi, below (/3 = 1.4, red lines) and above 
(/? = 1.6, black lines) the transition from diffusion to self- 
trapping. Here the initial state is a single-site (5-function with 
e = and A = 0.8. 



standing of the origin of the self-trapping in a lattice 
is based on energy conservation arguments Let us 
consider separately the contribution to the energy that 




FIG. 6: (color online) Critical value of the interaction strength 
for the transition to self-trapping, /3c, as a function of the dis- 
order strength. A, for diffusion from a single lattice site and 
for 6 = (blue squares) and = n (red circles). The red and 
blue dashed lines are the corresponding upper bounds for /3c 
obtained by calculating the bandwidth of the single-particle 
spectrum and using energy conservation arguments. The di- 
agram is schematically divided in three regions I, II and III. 
All states in I are self-trapped; in II, one finds diffusion, with 
soliton-like structures and discrete breathers; in III, the tran- 
sition from diffusive states to self-trapping strongly depends 
on the value of the phase 8 (i.e., the position of the initial 
wavepacket). The semi-axis A > 2 and (3 = corresponds to 
the regime of disorder induced localization for noninteracting 
particles. 



comes from the kinetic and potential terms together and 
the contribution that comes from the interacting term, 
H{t) = Ho{t) + Hint{t)- If the gas is subject to a pe- 
riodic potential in the tight-binding approximation and 
its dynamics is restricted to the lowest Bloch band, as 
supposed in deriving Eq. 12]), the term Hq in the Hamil- 
tonian is upper bounded. Let us call this upper bound 
j^max ^ ■\Yhenever the initial energy of the interacting sys- 
tem is larger than this upper bound, H{t = 0) > i?™"^, 
one can prove that the system cannot reach a situation 
where Hintit) — 0, at any t > 0, without violating the 
energy conservation. This means that, under these con- 
ditions, part of the interaction energy must be trapped 
in the system in the form of a localized peak that does 
not spread. In other words, whenever H{t = 0) > 
an initially localized wavepacket cannot spread to zero in 
the whole space. This argument, in general, does not pro- 
vide a precise quantitative estimate of the critical value 
(3c, but it gives a reasonable upper bound. 

Self-trapping of particles in a ID quasi-periodic poten- 
tial for A < 2 has been already discussed in R ef. |33l| (for 
the case of purely random potential see [tI. [Tol.lll|). Here 
we provide a more systematic calculation of (3c and we 
compare the diffusion from a single-site to the one from 
a Gaussian wavepacket. 

A signature of the presence of self-trapping is a satu- 
ration of the participation number P{t) that, for /3 > /3c, 
reaches an asymptotic finite value, due to the trapping 
mechanism occurring at the center of the wavepacket, 
while the width w{t) keeps increasing owing to the ex- 
panding tails [131 • An example of self-trapping transition 
is shown in Fig. [5l where we show the results obtained 
by solving Eq. ^ for diffusion from a single-site in a 
quasi-periodic potential with a = (\/5 — l)/2- In the 
figure one can see the typical change of behaviour that 
occurs when (3 crosses the critical value (3c. The same 
figure shows also the difference in the density distribu- 
tions at i = 1000: the central peak is well visible in the 
self-trapped state and it is absent in the diffusive state, 
while the lateral, low density tails are similar. 

By systematically looking at the numerical results for 
w{t), P{t), IV'jXOP the /? vs. A plane, we obtain the 
diagram shown in Fig. [6l The values of pc are repre- 
sented by red circles and blue squares for 9 = n and 0, 
respectively. In region I, above the red circles, all points 
correspond to self-trapped states. For A < 2 we find that 
the value of /3c is practically independent of the phase 9 
and decreases as A is increased. In region II, we observe 
diffusion, often accompanied by solitonic structures and 
discrete breathers eventually spreading. Similar struc- 
tures in the numerical solutions of Eq. (j4]) , for diffusion 
from a single-site and for A = 0, have been already found 
in Ref. [33|. For A > 2 we find that /3c is strongly 9- 
dependent. In the figure we show the results for the two 
limiting values 9 = and = tt; in particular, in re- 
gion III, we find that all states are self-trapped for 6* = 
while they are diffusive for 9 = tt. The semi-axis A > 2 
and /3 = corresponds to the regime of disorder induced 
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localization for noninteracting particles. 

The phase dependence of f3c for A > 2 can be quali- 
tatively explained by the energy conservation arguments 
already mentioned above. In particular, we numerically 
calculate the maximum energy, , in the lowest Bloch 
band of the noninteracting single-particle spectrum and 
we compare this value to the initial energy of the interact- 
ing system, which is given by H{t = 0) = A coa{6) + f3/2. 
The upper bound for the transition to self-trapping is 
then given by the condition H{t — 0) = i?™"^, which 
implies 

f3 = 2(£;^°^ ± A) 

where the plus and the minus signs holds for 6 = n and 
9 = 0, respectively. These two upper bounds are repre- 
sented by the blue and red dashed lines in Fig. [S) 

Fig.[6]shows that, in the case of diffusion from a single- 
site, the self-trapping mechanism plays a rather impor- 
tant and nontrivial role, leaving almost no space to the 
observability of the interplay between disorder and inter- 
action. The region were this interplay might be observed, 
namely for A > 2 and small /3, where one expects to see 
the destruction of localization due to interaction, it is 
also the region where the dependence on the phase 9 is 
the largest. Unfortunately, in typical experimental situa- 
tions with Bose-Einstein condensates, the phase 9 is not 
controllable. Moreover, in the experiments the initial dis- 
tribution of atoms in the lattice sites is more similar to a 
Gaussian than a (5-function. This suggests that, while the 
single-site diffusion is conceptually important and widely 
investigated from the theoretical viewpoint, the diffusion 
from a Gaussian is also interesting and worth exploring. 

By repeating the same calculation of /3c as before, but 
starting from a Gaussian wavepacket of width a, we find 
two main results: i) if a is of the order of 5 or more, the 
time evolution of the width, the participation number 
and the density distribution becomes almost independent 
of the phase 9 for all values of A and /?; ii) self-trapping 
is strongly suppressed, especially in the region A > 2. 
An example is shown in Fig. [7] for A = 2.5 and an initial 
Gaussian of width cr = 5. As one can see, self-trapping 
does not occur even for values of the interaction parame- 
ters of the order of /3 ^ 50. This fact is important in view 
of the discussion about the delocalization induced by the 
interaction, which is the subject of the next section. 

B. Delocalization induced by the interaction 

Let us now investigate the interplay of interaction and 
localization in the diffusion of a Gaussian wavepacket, in 
the regime where the noninteracting system is localized 
(A > 2) and for a choice of parameters that rules out 
self-trapping. In this case we find that, as the interac- 
tion is turned on, a wavepacket that was localized for 
P = starts to expand sub-diffusively. We observe an 
asymptotic growth of both the width w{t) and the par- 
ticipation number P(t), according to the following laws: 



w{t) ~ and P{t) - with 71,2 in the range (0, 0.5). 
The same effect has previously been reported for single- 
site diffusion in a purely random systems [1, [13, El • In 
the absence of self-trapping we find that the coefficients 
7i and 72 are nearly equal, therefore in the following we 
will use 7 w 71 w 72. An example is shown in Fig [7] 
where one sees the effects of the interaction on the diffu- 
sion of an initial Gaussian wavepacket with cr = 5, and a 
disorder strength just above the localization transition, 
A = 2.5. The noninteracting case, which remains local- 
ized, is compared with three different values of the inter- 
action parameter, /3 = 1, /? = 10 and /3 = 50. Already for 
(3=1 there is an evident delocalization and this effect 
increases as P is increased in the sense that 7 becomes 
larger and the delocalization takes place earlier. A very 
similar behaviour is obtained also for the diffusion from 
a single-site, provided the phase 9 and the interaction /3 
are chosen in such a way to avoid self trapping (e.g., in 
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FIG. 7: (color online) Time evolution of the width of the 
wavepacket w(t) (a) and of the participation number (b) for 
A = 2.5 and for an initial Gaussian wavepacket with cr = 5. 
We compare the noninteracting case, /3 = 0, with three inter- 
acting cases: /3 = 1, 10, 50. The black straight lines represent 
a guide to the eye. Their slope is 0.2, 0.3 and 0.34 and is the 
same in (a) and (b). These lines suggest that 71 ~ 72 and 
that 7 depends on (3 (see text) . 
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which is for /? = 10, we consider Gaussian wavepack- 
ets with cr = 2, 5, 6, 10, 14. We find that changing the 
width of the initial wavepacket does not affect the spread- 
ing mechanism, even in the Umit a — *■ of a (5-function 
wavepacket. In fact, one can see from this figure that 
there is no visible dependence of the asymptotic slope, 7, 
on the initial width. 



Time 



V. CONCLUSIONS 



FIG. 8: (color online) Time evolution of the width of a Gaus- 
sian wavepacket with ct = 5, for (3 = IQ and different values 
of the disorder strength A. 



region III of Figl6]with 9 ^ n). 

When the disorder strength A is increased the local- 
ization gets more robust, in the sense that the onset of 
sub-diffusive spreading takes place for later times and 7 
becomes smaller. For large A we reach a situation where 
the delocalization process is no longer observable within 
our simulation time. This is shown in Fig|51 where we 
compare the time evolution of a Gaussian wavepacket for 
fixed P and for increasing values of the disorder strength 
A. Finally, in Fig[5] we show the behaviour of w{t) as 
function of time for different values of the initial width 
of the wavepacket. In the left plot, which is for /3 — 2, 
we compare a 5-function wavepacket with some Gaus- 
sian cases with cr — 2, 6, 10, 14 while in the right plot. 
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FIG. 9: (color online) Width w{t) as a function of time for 
A = 2.5 and two different values of the interaction strength, 
(3 = 2 (left column) and /3 = 10 (right column). In the left 
plot we compare the time evolution for a single-site (5-function 
(red line) and for Gaussian wavepackets with different values 
of the initial width, a = 2 (green line), 6 (blue line), 10 (purple 
line), 14 (azure line) while in the right plot we consider only 
Gaussian wavepackets with a = 2 (green line), 5 (black line), 
6 (blue line), 10 (purple line), 14 (azure line). The curves 
exhibit the same asymptotic slope, which is 0.23 (left plot) 
and 0.3 (right plot). 



The role of inter-particle interaction in the destruction 
of Anderson localization is rather nontrivial and is the 
subject of current investigations both theoretically and 
experimentally. Here we have considered the case of the 
exponential localization of a wavepacket in a ID quasi- 
periodic potential. Our main goal was to fill the gap 
between what is predicted by the Aubry-Andre model 
for noninteracting particles in such a potential and what 
can be actually observed in realistic experiments with 
ultracold bosons in bichromatic optical lattices. As a 
first step in this direction, we have studied the diffusion 
of a noninteracting wavepackets in a commensurate (pe- 
riodic) lattice and we have compared it with the case 
of an incommensurate (quasi-periodic) lattice. We have 
shown that the spatial periodicity of the commensurate 
lattice plays a key role in determining the type of ap- 
proach to the quasi-periodic limit in a sequence of com- 
mensurate approximants. This part of our analysis con- 
firms that the transition from diffusion to localization 
observed in Ref. [2] can correctly be interpreted in terms 
of the predictions of the Aubry-Andre model. A sec- 
ond step consists of including the atom-atom interaction. 
To this purpose we have numerically solved a discrete 
non-linear Schrodinger equation, which generalizes the 
Aubry-Andre model by introducing the interaction at the 
mean-field level. We have simulated the dynamics of mat- 
ter waves starting from either a ^-function localized in a 
single lattice site or a Gaussian wavepacket. In the for- 
mer case, we have found that the dynamics is dominated 
by self-trapping processes in a wide range of parameters, 
even for weak interaction. Conversely, in the latter case, 
self-trapping is significantly suppressed and the destruc- 
tion of localization by interaction is more easily observ- 
able. In particular, we find that Gaussian wavepackets, 
which remain localized for noninteracting particles, start 
to spread sub-diffusively in the presence of interaction. 
This result is consistent with previous predictions for in- 
teracting particles in purely random potentials. We have 
systematically investigated the transition from localiza- 
tion to diffusion as a function of the strength of both the 
interaction and the disorder. Our analysis is intended to 
stimulate further experimental work on the diffusion of 
atomic matter waves in bichromatic lattices. 
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